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ABSTRACT 


Solar energetic particles (SEPs) affect the solar-terrestrial space environment 
and become a very important aspect in space weather research. In this work, 
we numerically investigate the transport processes of SEPs in three-dimensional 
interplanetary magnetic held, with an emphasis on the longitudinal distribution 
of SEPs in the heliosphere. We conhrm our previous Ending that there exists an 
east-west longitudinal asymmetry in the SEP intensities, i.e., with the same lon¬ 
gitude separations between the solar source centers and the magnetic footpoint 
of the observer, the huxes of SEP events originating from solar sources located on 
the eastern side of the nominal magnetic footpoint of the observer are systemati¬ 
cally larger than those of the SEP events originating from sources located on the 
western side. We discuss the formation mechanism of this phenomenon, and con¬ 
clude that the longitudinally asymmetric distribution of SEPs results from the 
east-west azimuthal asymmetry in the topology of the heliospheric magnetic held 
as well as the ehects of perpendicular dihusion on the transport of SEPs in the 
heliosphere. Our results will be valuable to understanding Sun-Earth relations 
and useful for space weather forecasting. 

Subject headings: interplanetary medium - magnetic fields - solar-terrestrial relations 
- Sun: flares - Sun: particle emission 
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1. Introduction 


Solar energetic particles (SEPs), which are charged energetic particles occasionally 
emitted by the Snn, risk the health of astronants working in space and damage electronic 
components on satellites, so they have become a very important aspect affecting solar- 
terrestrial space environment and space weather. Theoretically, SEPs observed in the 
interplanetary magnetic held provide fnndamental information regarding the acceleration 
mechanisms and transport processes of charged energetic particles. Therefore, the snbject 
has become a focns of astrophysics, space physics, and plasma physics. 


Throngh several decades of investigations in the commnnity with spacecraft 
observations and theoretical modeling, signihcant progresses have been achieved toward 
a bette r nnde rstanding of the transport processes of SEPs. As a pioneering work. 


Parker! fllQGSll provided a diffnsion eqnation to investigate the modnlation of galactic 


cosmic rays and the transport of SEPs. It is very difhcnlt to analytically solve the 


adopted fe.g.. 

Ne; & Wone 

1979; 

Rnholo 

1991; 

Zhane; et al. 

2009; 

He & Oin 

2011; 

He et al.l 

2011; 

Giacalone & Jokiph 

2012) 

. Jokipi] 

( 

196( 

3) proposed the well-known and classical 


qnasi-linear theory (QLT) to describe the diffnsion of charged particles in a tnrbnlent 
magnetic held. The theoretical investigations and observational analyses based on SEP 


angle in the particle pitch 


angle distri 


Dntions 


Beeck & Wibberenz 

1986; 

Beeck et al. 

1987) 




Hasse 

mann & Wibberenz 

1968, 

1970 


Schlickeiserl fl2002h laid complete physical 


fonndatio ns and presentee 


research. 


Matthaens et al 


qrofqnnd mathematical techniqnes for cosmic ray transport 


(120031 ) provided a nonlinear gniding center (NLGC) theory for 


charged particles’ 


of test particles. 


perpendicnlar dihu sion to acconnt for the nnmerical simulation resnlts 


Shalchi et al. 


(1200411 presented analytical forms for the resnlts from the 
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NLGC theory. IZank et al.l fj2000l) and lLi et al.l (120031) provided a dynamical time-dependent 
model of charged particle acceleration at a propagating and evolving interplanetary shock. 


Recently, IZhang et al.l (120091) presented a model calcnlation of SEP propagation in 


a three-c 


He et ah 


i mens ional interplanetary magnetic held with the effect of magnetic tnrbnlence. 


(120111) investigated the effects of particle sonrce characteristics on SEP observations 


at 1 AU and fonnd that the perpendicnlar diffnsion has a very important inhnence on 
the propagation of SEPs in the heliosphere, particnlarly when a spacecraft is not directly 
connected to the acceleration regions either on the Snn or near the CME-driven shocks 
by the interplanetary magnetic held lines; in snch cases, the earliest arriving particles 


He & Wan 

(2012a 

1 presented a 

spatially dependent para! 

el and 

Furthermore, 

He & Wan 

(2ni2bll 


2012bl) provided a direct approach to explicitly investigate the 


spatially dependent perpendicnlar mean free path of SEPs in a tnrbnlent and spatially 
varying magnetic held. They reported that for physical conditions representative of the 
solar wind and the magnetic tnrbnlence, the ratio Ax/A|| of the perpendicnlar to the parallel 
mean free path remains in the range 0.01 — 0.20; however, when the tnrbnlence strength 
5B/B is snfhciently large, the ratio A_l/A|| wonld approach or exceed nnity. 


Observationally, 


Cane et al. 


(1198811 investigated ti me proh l es of 2 35 intense proton 


events to docnment the longitnde-dependent prohles. 


ReamesI (119991) presented some 


observational evidences of longitndinal ehects on the transport processes of SEPs. They 
fonnd that the interplanetary spiral magnetic held causes an asymmetry in the intensity 
time prohles of SEP events originating from eastern and western longitudes on the solar 
disk. As we know, SEP s are usually p roduced by solar hares or/and r elated to coronal mass 


ejections (CMEs) (e.g.. 


Reames 


I 999 I) . Recently, 


Dresing et al. 


(I2012i) investigated in detail 
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the large longitudinal spread of SEPs during the 17 January 2010 solar event and suggested 
with observational evidence that the wide longitudinal spread of particles is mainly due to 
perpendicular diffusion in the interplanetary medium. 


In 


He et ah 


(1201 ih . the authors found that with the same separation in heliographic 


longitude between the magnetic footpoint of the observer and the solar sources, the SEPs 
associated with the sources located east are detected earlier with larger fluxes than those 


associated with the sources located west. Th is interest i ng p 


east-west azimutha l asymmetry of SEPs” by 


He et al. 


lenomenon was called the 


(120111 ). The simulation results in 


Giacalone fc Jokipiil (120121) also showed similar asymmetric features of SEP distribution. 


In this paper, we focus on the so-called east-west azimuthal asymmetry of SEPs in view 
of numerical investigations. There a 
dependence of SEP intensities (e.g.. 


Cane et al. 


Reames 

1999 ) 


19991) . Basically, however. 


these previous studies only analyzed the effects of heliographic longitude of particle sources 
on the intensity profiles of SEP events. In our simulations, we mainly investigate the 
relatively longitudinal distribution of SEP intensities in the heliosphere, and find the 
longitudinally asymmetric distribution phenomenon. With the same longitude separations 
between the solar source centers and the magnetic held line footpoint of the observer, the 
huxes of SEP events originating from solar sources located on the eastern side of the nominal 
magnetic footpoint of observer are systematically larger than those of the SEP events 
originating from sources located on the western side. We conclude that the longitudinally 
asymmetric distribution of SEPs results from the east-west azimuthal asymmetry in the 
geometry of the heliospheric magnetic held as well as the ehects of perpendicular dihusion 
on the transport processes of SEPs in the heliosphere. Our results will be valuable for 
understanding solar-terrestrial relations and useful in space weather forecasting. Some 
investigation results related with the statistical and numerical studies on the longitudinally 
asymmetric distribution of SEPs were reported in the 33rd International Cosmic Ray 
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Conference flHe &: Wan 


2ni8ah 


This paper is organized as follows. Firstly, in Section 2, we present a five-dimensional 
focussed transport model containing most of the usually used transport mechanisms. The 
simulation method for numerically solving the focussed transport equation is also provided 
in this section. In Section 3, the simulation results of the relatively longitudinal distribution 
of SEPs will be presented. We present some implications of our simulation results in Section 
4 and a discussion on the formation mechanism in Section 5. A summary of our results will 
be provided in Section 6. 


Numerical Model and Method 


In our model, the five-dimensional focussed transport equation that governs the 

of SEPs can be written as (e.g.. 


Schlickeiser 

2002 ; 

Zhang et al. 

2009 


He et al 


2011 


Qin et al. 


201 ih 


df ^ df 

-A A A 

dfi \ d/i 


dt dp dt dp 


dx dx ) 




( 1 ) 


where x is particle’s position, z is the coordinate along the magnetic field spiral, p is 
particle’s momentum, p is the pitch-angle cosine of particle, t is time, v is the particle 
velocity, is the solar wind speed, and Q is the source term. In the transport model, 
we directly input a source of accelerated particles as a product of either a solar flare or 
the CME-driven shock. We focus on the high-energy SEPs accelerated and released into 
the system near the Sun. The particle injection can be seen as a point source in the radial 
direction. For observations far out near 1 AU or larger radial distances, this assumption is 


parti cularly valid and is commonly used in the study of SEP transport (e.g.. 


Zhang et al 


20091 ). The particle momentum change term dp/dt in Equation ([T]) due to the adiabatic 
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cooling effect can be written as fjSkillindll97lll) 


dp 

Tt=-^ 


1 - fdV^ 


+ 


y 


+ p 


>dVS 


2 \ dx dy J dz 

In addition, dp/dt, which includes the magnetic focusing effect and t 


( 2 ) 


solar wind flows effect, can be expressed as (e.g., 
19971 1 


Roelof 


19691: 


le divergence of the 


Isenberg 


1997 : 


Kota fc Jokipii 


dp 

dt 


2 r 


1 -h 
2 

l-p2 r 


V OB 
'B~d^ 

V 


+ h 
aKf 


X 

dx 


+ 


gys^ 

dz 


dy 


+ 




dx ' dy ~ dz J 
where B is the background interplanetary magnetic held with direction z, and the magnetic 
focusing length L is dehned as L = (z ■ yin5)“^. Therefore, the adiabatic cooling effect 
and the divergence of the solar wind hows ehect generated by the solar wind velocity 
components, expressed as Equations ([2]) and ([3]), respectively, have been included in our 
model. The interplanetary magnetic held in the inner heliosphere is too strong and the 


particle energies in typical SEP events are not 


ligh enough, so in this work we neglect the 


particle drift in the nonuniform magnetic held (IZhang et ah 


20091). 


Under the dihusion approximation for a nearly isotropic pitch-angle distribution, the 


Earl 

1974; 

He & Schlickeiser 

2Q14) 


hlokipii 

1999; 

Hasselmann & Wibberenz 

1998 


A„ = ^ 


r- + l 


2\2 


(1-h^) 


-dp. 

Accordingly, the radial mean free path can be dehned as 

Xr = A|| cos^ ^|J, 

where 'd is the angle between the local magnetic held direc 


Equation ([5]) , cos^ can be written as flNg fc Gleeson 


direct 

ion and the radia 

1971; 

He & Wan 

2012 ai 


( 4 ) 

( 5 ) 

direction. In 


cos^ ^ sin^ 9 ) , 


( 6 ) 













































































where is the solar wind speed, hi is the angular rotation velocity of the Sun, and r and 
6 are the coordinates of the heliocentric spherical coordinate system {r,6,(j)), namely, r is 
the heliocentric radial distance, and 6 is the colatitude, which is measured from the rotation 
axis of the Sun. 


We use a form of the pitch-angle diffusion coefficient as fe.g.. lBeeck fc Wibberen7j|l986l: 


He et af 


201 ih 


Dl, = DJ ca? i< = DovR-^/^ (|f,r' + ft) (1 - P), 


( 7 ) 


where Dq is a constant indicating the magnetic turbulence strength and R is the particle 
rigidity. The constant h is needed to simulate the particles’ ability to scatter through /i = 0. 
In order to simulate the nonlinear effect to cause large at /i = 0, we use a relatively 
large value of h = 0.2. The constant q is related to the power spectrum of the magnetic field 
turbulence in the inertial range, chosen to be 5/3 in the model. Additionally, we assume 
that the two perpendicular diffusion coefficients, A^, and A^, are the same and independent 


of yU. According to the observational de terminations and 


mean free paths in previous works (see 


He fc Wan 


heoretical investigations of the 


2013bl. and references therein), we 


typically use the radial mean free path A^ = 0.3 AU (corresponding to the parallel mean 
free path Ay = 0.67 AU at 1 AU) and the perpendicular mean free paths \x = \ = 0.006 
AU. Actually, altering the values of the parameters in the model will not qualitatively 
change the simulation results, i.e., the longitudinally asymmetric distribution of the SEP 
fluxes, which will be discussed later. 


In order to numerically solve the focused transpor t equation 


time-backward Markov stochastic process method (e.g.. 


Zhang 


19991) 


we utilize the 
This approach could 


deal conveniently with an expanded source energy spectrum. When using the time-backward 
stochastic process method, we trace SEPs back to the initial time, and only those particles 
in the source region at the initial time contribute to the statistics. The five time-backward 
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stochastic differential equations transformed from the focused transport equation ([T]) can 
be described as follows: 


dX = V^dW^is)-V^^^ds 
dY = .y^dWyis)-Vy^^ds 
dZ = -{fiV + V^^^)ds 


d/i = ^/2Df,ydWf,{s) 


1 - h' 
2 

dD 


2 r 


y 


MM 


+ 


2D 


1111 


9/i M + /i 


dV^'^ 

dx 

ds 


+ 


svr _ /r. 


dy 


dz 


ds 


dP = P 




dVr ^ dVy^" 


+ jj 


,dvr 


ds, 


( 8 ) 


dx dy / dz 

where (X, Y, Z) is the pseudo-position, V is the pseudo-velocity, P is the pseudo¬ 
momentum, and lyE(^), iy 2 ;(^), and are Wiener processes. In our simulations, the 

Parker interplanetary magnetic field B is set so that its magnitude at lAU is 5nT, and the 
solar wind speed is typically set as = 400 km s~^. 


The source term Q in Equation ([T]) with a power-law spectrum 7 is assumed as flReid 


196411 


Q{z < O.OSkiy Ek, 9, 0, t) = y exp (~ - —^ ^{9, 0), (9) 

t \ t TlJ 

where ^{9,(j)) is a function of the heliographic latitude 9 and longitude 0 that describes the 
spatial distribution of SEP source strength, Ek is the kinetic energy of the source particles, 
and Tc and tl denote the rise and decay timescales of source release profile, respectively. 
Although the injection form in Equation ([9]) is originally based on an assumption of lateral 
diffusion in the corona, it can be used as a short pulse of source particle injection, especially 
when the shock is near the Sun. Typically, the values of Tc and tl are a few hours, which 
are much sh orter than the tirn escale to reach the maximum flux or the decay phase of 


SEP events (IZhang et ah 


20091 ). Therefore, Tc and tl only affect the flux profiles at the 
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beginning, but not the long-term when t S> tl- In this work, we typically set ■y = 3, Tc = 0.1 
days and tl = 0.25 days. As we know, typical sizes of SEP s ources (flares or CMEs) are tens 


of deCTees wide in helio graphic latitude and longitude (e.g.. 


2001 


Hundhausen 


Wang et ah 


1993; 


Maia et al. 


20061) . It is not easy to hnd out the exact sizes of SEP sources; in the 


simulations, we set SEP sources at r = 0.05A[/ with limited coverage of 30° in latitude and 
longitude. Furthermore, we assume a uniform spatial distribution in the source regions. In 
addition, an outer absorptive boundary of pseudoparticles is set at r = 50A17. 


To obtain the SEP flux for each case, we simulate 6 x 10^ pseudoparticles on a super¬ 
computer cluster with message passing interface. Generally, the unit of omnidirectional flux 
is used as — sr~^ — — MeV~^] in this work, however, we use an arbitrary unit for 

convenience in plotting hgures. 


3. Numerical Results 

For each particle energy investigated, we simulate 18 cases corresponding to different 
source locations to show the effect of solar source location on the SEP flux observed in 
the interplanetary space. All the solar particle sources have the same coverage, i.e., 30° in 
latitude and longitude, but their centers are located at various heliographic longitudes in 
the solar equator or at 90° colatitude. Specihcally, the longitude separations between the 
centers of the 18 solar sources and the magnetic footpoint of the observer are set to be the 
following values: 0°, 20° west, 20° east, 40° west, 40° east, ..., 160° west, 160° east, 180° 
east. All the other conditions and parameters, such as the source particle intensity, solar 
wind speed, turbulent magnetic held, etc., are the same for each case. In the simulations, 
we typically investigate 20 MeV, 50 MeV, and 100 MeV protons detected at 1 AU, 90° 
colatitude. 
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A diagram, i.e., Figure [H serving as an example, is presented to illustrate the 
longitudinal locations of the particle sources on the solar surface relative to the magnetic 
footpoint of the observer, which is located at 90° colatitude. In the scenario of the diagram, 
we can see that the center of the middle source on the Sun is connected directly to the 
spacecraft by interplanetary magnetic held line, whereas the centers of the solar sources on 
the left and right are 60° east and 60° west, respectively, away from the magnetic footpoint 
of the observer. 

Figure 12] shows the hux-time prohles of the 20 MeV SEPs originating from the different 
solar sources with different longitudinal locations relative to the magnetic footpoint of 
the observer at 1 AU, 90° colatitude. In Figure IH the solid and dashed lines denote the 
SEP omnidirectional huxes originating from the east and west solar sources, respectively. 
The different colors of the curves represent different longitudinal distances between the 
magnetic footpoint of the observer and the solar sources. The 0° relatively longitudinal 
distance denotes the case where the nominal magnetic footpoint of the observer and the 
solar source center locate at the same heliographic longitude. As we can see, the SEP 
flux in this case is the largest among all the cases investigated. This means that when 
the observer in the interplanetary space is nearly directly connected to the solar source 
by heliospheric magnetic held lines, the SEP hux is larger than that otherwise. However, 
a careful investigation based on our numerical simulations indicates that the maximum 
SEP hux shifts from the location of best magnetic connection and toward the east SEP 
sources with centers on relative longitudes ~ —0.7°— ~ —1.5°. We note that the shift 
extents depend on the ratio A_l/A|| of the perpendicular to the parallel mean free paths, the 
geometry of the magnetic spiral, the energy of SEPs, the coverage of the SEP sources, and 
the particle spatial distribution in the sources. In Figured we can also see that the farther 
the magnetic held footpoint of the observer is away from the solar source, the smaller is the 
SEP hux observed and also the later the SEP event arrives at the observer’s position (see 
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also 


He et al 


20111). Moreover, with the same longitude separations between the magnetic 


held footpoint of the observer and solar sources, the SEP huxes from solar sources located 
on the eastern side of the observer footpoint are larger than those from solar sources 
located on the western side, and the times of onsets and maximum huxes of SEP events 
from east sources are earlier than those from west sources. In other words, there exists an 
east-west longitudinally asymmetric distribution of SEPs in the heliosphere. In addition, 
we can obviously observe that in the decay phase the hux prohles of all the SEP cases have 
almost the same decay rate, which is called the SEP reservoir phenomenon in spacecraft 
observations. Therefore, our numerical simulations with a series of diherent SEP cases 
almost covering 360° in heliographic longitude successfully reproduce the so-called SEP 
reservoir. The SEP reservoirs are observed by spacecraft at both low and high heliolatitudes, 
and by spacecraft at very different heliolongitudes and radial distances. In our simulations 
with perpendicular diffusion, we reproduce various SEP reservoirs (some results not shown 
here) at different heliospheric locations (longitude, latitude, and radial distance) without 
invoking the hypothesis of a reflecting boundary or diffusion barrier. Actually, it is difficult 
to imagine such an “overwhelming” reflecting boundary or diffusion barrier covering all the 
longitudes, all the latitudes, and even all the radial distances. Therefore, the perpendicular 
diffusion should be the dominant, if not the only, physical mechanism responsible for the 
formation of the SEP reservoirs. 


The circles with different colors on the flux profiles in Figure |2] denote the peak fluxes 
of the corresponding SEP cases. We extract the information of the peak fluxes and the 
relevant longitudinal distances in the SEP cases investigated from Figure Ej and present 
it separately. In Figure [3l we show the peak fluxes of the 20 MeV SEPs originating from 
the different solar sources with different longitudinal locations relative to the magnetic 
footpoint of the observer at 1 AU, 90° colatitude. The X-axis in Figure |3] is the relatively 
longitudinal distance between the solar source center and the magnetic footpoint of 
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the observer. The vertical dashed line in the middle of Figure [3] denotes the 0° relatively 
longitudinal distance, where the nominal magnetic footpoint of the observer and the solar 
source center locate at the same heliographic longitude. As we can see, the peak flux in 
the SEP case with about 0° relatively longitudinal distance is the largest among the cases 
investigated. This means that when the observer in the interplanetary space is nearly 
directly connected to the solar source by heliospheric magnetic field lines, the SEP peak flux 
is larger than that otherwise. However, a careful investigation indicates that the maximum 
SEP peak flux shifts from the place of best connection and toward the east SEP sources 
with centers on relative longitudes ~ —0.7°— ~ —1.5°. We connect the peak fluxes of SEPs 
in the simulated cases on the left-hand and right-hand halves with the red and blue curves, 
respectively. It can obviously be seen from Figure |3] that the farther the solar source center 
is away from the magnetic footpoint of the observer, the smaller is the peak flux of SEPs 
observed. We further mirror the peak fluxes of SEPs on the right-hand half to the left-hand 
half with a blue dotted curve. We can obviously observe that with the same longitude 
separations between the solar source centers and the magnetic footpoint of the observer, the 
peak fluxes of SEP events originating from solar sources located on the eastern side of the 
observer footpoint are systematically larger than those of the SEP events originating from 
sources located on the western side. Therefore, our relatively complete model calculation 
of SEP propagation with the effect of perpendicular diffusion shows that there exists an 
east-west asymmetry in the SEP intensities. 

Figure m shows the peak fluxes of the 50 MeV SEPs originating from the different 
solar sources with different longitudinal distances relative to the magnetic footpoint of the 
observer at 1 AU, 90° colatitude. As in Figure |3l we can also clearly observe that with the 
same longitudinal distances between the solar source centers and the magnetic footpoint 
of the observer, the peak fluxes of the SEP events originating from solar sources located 
on the eastern side of the observer footpoint are systematically larger than those of the 
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SEP events originating from sources located on the western side. Therefore, the numerical 
simulation of the transport process of 50 MeV SEPs also shows that there exists an east-west 
longitudinally asymmetric distribution of SEPs in the heliosphere. We note that the subtle 
difference between Figure [3] and Figure 0] results from the interplay of diffusion coefficients 
K± and K\\ as well as the weak statistics of SEP fluxes in the cases with long longitudinal 
distances between the solar sources and the magnetic footpoint of the observer. Figure 0] 
shows the peak fluxes of the 100 MeV SEPs originating from the different solar sources 
with different longitudinal distances relative to the magnetic footpoint of the observer at 1 
AU, 90° colatitude. It can also evidently be seen that with the same longitudinal distances 
between the solar source centers and the magnetic footpoint of the observer, the peak fluxes 
of the SEP events originating from solar sources located on the eastern side of the observer 
footpoint are systematically larger than those of the SEP events originating from sources 
located on the western side. 


4. Implications of Simulation Results 


Generally, an SEP event with a higher peak flux will reveal a larger intensity relative to 


other events through the entire evolution process (see Figure [2] and the results in 


He et ah 


(1201 ih ). In general, the probability of the SEP events being observed by spacecraft near 
Earth’s orbit depends on the intensities of the particles in the events. An intense SEP event 
with high flux will experience relatively weak influences from interplanetary structures (such 
as magnetic cloud, heliospheric current sheet, corotating interaction region, etc.) during the 
propagation process in the heliosphere. However, the not so intense SEP events with small 
fluxes will be significantly affected by the interplanetary structures and then to weaken 
or even dissipate before entering Earth’s orbit. For example, numerous observations have 
shown that a magnetic cloud can considerably affect the propagation conditions of SEPs 




15 


and low-energy cosmic r ays, which is 


energetic particles (e.g., 


Blanco et al 


comm only observed as a decrease in the flnxes of the 


20131) . Conseqnently, intense SEP events wonld have 


a larger probability of reaching Earth and being observed by spacecraft at 1 AU; whereas 
the relatively weak SEP events wonld have a lower probability to arrive at Earth and to be 
recorded by spacecraft. 


As shown in the simulation results above, the SEP events originating from solar 
sources located on the eastern side of the magnetic footpoint of the observer reveal larger 
fluxes than those originating from solar sources located on the western side. Therefore, 
the SEP events originating from solar sources located at eastern relative longitudes would 
have a higher probability to reach Earth and to be observed by spacecraft near 1 AU. On 
the contrary, the SEP events from solar sources on the west side relative to the magnetic 
footpoint of observer will be more easily influenced by the interplanetary medium and 
structures, and then to weaken or even dissipate during their propagation processes. As 
a result, the spacecraft at 1 AU would have a larger probability to miss the SEP events 
originating from solar sources located at western relative longitudes. 


5. Discussion 


Generally, as viewed from above the north pole of the Sun, the Parker spiral magnetic 
held lines originating from the solar surface would meander clockwise to somewhere in the 
interplanetary space with more eastern longitudes than their footpoints on the solar surface. 
This scenario can be seen in Figure |6l As we can see, there is an east-west azimuthal 


asymmetry in the topology of the Parker interplanet ary magnetic he^ 


the theoretical and obse rvational investigations (e.g.. 


Matthaeus et al. 


d. A c cording to 


2003 


Bieber et al. 


2004 


He fc Wan 


2012 bl) . for physical conditions representative of the solar wind, the 


perpendicular diffusion coefficient would be a few percent of the parallel diffusion coefficient 
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in magnitude. In general, SEPs would mainly transport along the average interplanetary 
magnetic field after they were produced from the surface of the Sun or a CME-driven 
shock. Therefore, when the observer is nearly directly connected to the solar source by 
the interplanetary magnetic held lines, the omnidirectional hux and peak hux of SEPs are 
larger than those in the cases where the observer is not connected directly to the solar 
sources. In general, the farther the solar source center is away from the magnetic footpoint 
of the observer, the smaller are the omnidirectional hux and peak hux of SEPs observed. In 
addition to the parallel dihusion, the SEPs will experience perpendicular dihusion to cross 
the heliospheric magnetic held lines during their transport processes in the interplanetary 
space. This is why the SEP events can be observed by a spacecraft that is not directly 
connected to the acceleration regions. In the limit of a point source release of SEPs and no 


perpendicular dihusion, the particles would a. 
held line rotating with the Sun. However, as 


so transport in long i tude a s a result of the 


Giacalone fc Jokipiil (120121) pointed out, in 


this case, the intensity prohle of the event observed near 1 AU would be a spike of SEPs 
only seen when the held line containing the particles crosses the observer and would be 
zero for all other times, which is inconsistent with the usually observed shape of intensity 
prohle of SEPs. A careful investigation based on our numerical simulations indicates that 
the maximum SEP hux shifts from the location of best magnetic connection and toward 
the east SEP sources with centers on relative longitudes ~ —0.7°— ~ —1.5°. Note that the 
shift extents depend on the ratio Ax/Ay, the geometry of the magnetic spiral, the energy of 
SEPs, the coverage of the SEP sources, and the particle spatial distribution in the sources. 


Due to the ehects of the perpendicular dihusion and the azimuthally asymmetric 
geometry of the heliospheric magnetic held, with the same longitudinal distances between 
the solar sources and the magnetic footpoint of the observer, the SEPs originating from 
sources located on the eastern side of the observer footpoint are easier and more frequent 
than those from sources located on the western side to arrive at the observer in the 
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interplanetary space. Consequently, with the same longitudinal distances between the solar 
sources and the observer footpoint, the fluxes of the SEP events associated with the solar 
sources located east are systematically larger than those of the SEP events associated with 
the sources located west. Accordingly, the peak intensities of the SEP events originating 
from solar sources located on the eastern side of the observer footpoint are also larger than 
those of the SEP events originating from sources located on the western side. Therefore, we 
propose that the longitudinally asymmetric distribution of SEPs results from the east-west 
azimuthal asymmetry in the topology of the Parker interplanetary magnetic held as well as 
the effects of perpendicular diffusion on the transport processes of SEPs in the heliosphere. 
Furthermore, the interplay between the perpendicular diffusion (A_l) and the parallel 
diffusion (A||) plays a very important role in determining the azimuthal asymmetry extents 
of the SEP distribution in the interplanetary space. 


6. Summary and Conclusions 


In the previous work of 


He et ahl (120111) . the authors reported and predicted that with 


the same longitude separation between the magnetic footpoint of observer and the solar 
sources, the SEPs produced and released from the sources located east are detected earlier 


lis interes 


He et ah 


IM. 


with larger fluxes than those associated with the sources located west. T 
phenomenon was called the “east-west az imuthal asymmetry of SEPs” by 
The figures in iGiacalone fc Jokipiil (120121) using Parker’s equation also displayed east-west 


( 120111 ) . 


asymmetry in the SEP distribution. In this work, we systematically study the longitudinally 
asymmetric distribution of SEPs in view of numerical investigations. We carry out a series 
of relatively complete model calculations of SEP propagation with perpendicular diffusion 
in a three-dimensional interplanetary magnetic field. Our simulation results show that with 
the same longitude separation between the solar sources and the magnetic footpoint of 
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the observer, the flux of SEPs released from the solar source located east is systematically 
higher than that of SEPs originating from the source located west. We conclude that 
the longitudinally asymmetric distribution of SEPs results from the east-west azimuthal 
asymmetry in the topology of the Parker interplanetary magnetic held as well as the effects 
of perpendicular diffusion on the transport processes of SEPs in the heliosphere. The 
interplay between the perpendicular diffusion (A^) and the parallel diffnsion (A||) plays 
an important role in determining the extents of the azimnthal asymmetry of the SEP 
distribntion in the heliosphere. 

In general, the probabilities of the SEP events being observed by spacecraft in 
the interplanetary space are proportional to the particle flux intensities in the events. 
Conseqnently, the SEP events originating from solar sonrces located at the eastern relative 
longitndes would be detected and recorded with somewhat higher probabilities by the 
observer near 1 AU. Therefore, the results presented in this work will be valuable for 
understanding the solar-terrestrial relations and usefnl in space weather forecasting. 
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Observer 


Fig. 1.— Diagram to illustrate the longitudinal locations of the particle sources in the 
solar equator relative to the magnetic footpoint of the observer, which is also located at 90° 
colatitude. The center of the middle solar source is connected directly to the spacecraft by 
the interplanetary magnetic held line, whereas the centers of the solar sources on the left 
and right are 60° east and 60° west, respectively, away from the magnetic footpoint of the 
observer. 



Omnidirectional Flux 
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Fig. 2.— Flux-time profiles of the 20 MeV SEPs originating from the different solar sources 
with different longitudinal locations relative to the magnetic footpoint of the observer at 
1 AU, 90° colatitude. With the same longitude separations between the magnetic held 
footpoint of the observer and solar sources, the SEP huxes from solar sources located on the 
eastern side of the observer footpoint are larger than those from solar sources located on the 
western side. In addition, the SEP reservoir phenomenon is reproduced. 













Flux maximum 
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Fig. 3.— Flux maxima (peak fluxes) of the 20 MeV SEPs originating from the different 
solar sources with different longitudinal locations relative to the magnetic footpoint of the 
observer at 1 AU, 90° colatitude. With the same longitude separations between the solar 
source centers and the magnetic footpoint of the observer, the flux maxima of the SEP events 
originating from solar sources located on the eastern (left) side of the observer footpoint are 
systematically larger than those of the SEP events originating from sources located on the 
western (right) side. 
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Fig. 4.— Flux maxima (peak fluxes) of the 50 MeV SEPs originating from the different 
solar sources with different longitudinal locations relative to the magnetic footpoint of the 
observer at 1 AU, 90° colatitude. With the same longitude separations between the solar 
source centers and the magnetic footpoint of the observer, the flux maxima of the SEP events 
originating from solar sources located on the eastern (left) side of the observer footpoint are 
systematically larger than those of the SEP events originating from sources located on the 
western (right) side. 
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Fig. 5.— Flux maxima (peak fluxes) of the 100 MeV SEPs originating from the different 
solar sources with different longitudinal locations relative to the magnetic footpoint of the 
observer at 1 AU, 90° colatitude. With the same longitude separations between the solar 
source centers and the magnetic footpoint of the observer, the flux maxima of the SEP events 
originating from solar sources located on the eastern (left) side of the observer footpoint are 
systematically larger than those of the SEP events originating from sources located on the 
western (right) side. 
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Fig. 6.— Sketch to show the east-west azimuthal asymmetry in the topology of the Parker 
interplanetary magnetic held. Due to the effects of the perpendicular diffusion and the 
azimuthally asymmetric geometry of the interplanetary magnetic held, with the same lon¬ 
gitudinal distances between the solar sources and the magnetic footpoint of the observer, 
the SEPs originating from sources located on the eastern side of the observer footpoint are 
easier and more frequent than those from sources located on the western side to arrive at 
the observer in the interplanetary space. 





